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ABSTRACT 

Context. Over the past decade, direct imaging has confirmed the existence of substellar companions on wide orbits from their parent 
stars. To understand the formation and evolution mechanisms of these companions, their individual, as well as the full population 
properties, must be characterized. 

Aims. We aim at detecting giant planet and/or brown dwarf companions around young, nearby, and dusty stars. Our goal is also to 
provide statistics on the population of giant planets at wide-orbits and discuss planet formation models. 

Methods. We report the results of a deep survey of 59 stars, members of young stellar associations. The observations were conducted 
with the ground-based adaptive optics system VLT/NaCo at L'-band (3.8/jm). We used angular differential imaging to reach optimal 
detection performances down to the the planetary mass regime. A statistical analysis of about 60 % of the young and southern A-F 
stars closer than 65 pc allows us to derive the fraction of giant planets on wide orbits. We use gravitational instability models and 
planet population synthesis models following the core-accretion scenario to discuss the occurrence of these companions. 
Results. We resolve and characterize new visual binaries and do not detect any new substellar companion. The survey's median 
detection performance reaches contrasts of 10 mag at 0.5" and 11.5 mag at 1.0". We find the occurrence of planets to be between 
10.8 and 24.8 % at 68 % confidence level assuming a uniform distribution of planets in the interval [1, 13] Mj and [1, 1000] AU. 
Considering the predictions of planetary formation models, we set important constraints on the occurrence of massive planets and 
brown dwarf companions that would have formed by gravitational instability. We show that this mechanism favors the formation 
of rather massive clump {M^t^^p > 30 Mj) at wide (a > 40 AU) orbits which might evolve dynamically and/or fragment. For the 
population of close-in giant planets that would have formed by core accretion (without considering any planet - planet scattering), 
our survey marginally explore physical separations (< 20 AU) and cannot constrain this population. We will have to wait for the next 
generation of planet finders to start exploring that population and even for the extremely large telescopes for a more complete overlap 
with other planet hunting techniques. 

Key words, instrumentation : adaptive optics - stars : young, nearby, dusty - methods : statistical - planetary system 

1 . Introduction iMavor et all (l201lh find that 50 % of solar-type stars harbor at 

least one planet of any mass and with period up to 100 days. This 

Most of the giant planets have been discovered so far thanks to occuiTence decreases to 14 % when considering giant planets 

indirect techniques (radial velocity and transit) at short orbits (< j^^gg^ j^an 0.3 Mjand vari es if we consider giant planets around 

5 AU). Almost 20 years of systematic seai'ch lead to numerous lower/higher mass stars dCumminge^ Si LTohnson et al.l 

surveys around solar-type, lower/higher mass (|Endl et al. || 2006|: 2OIO; Mayor et al. mil) (see Table[D. These rates thus confirm 

iBonfils et alj|20 131; 'Lagrange et al. 2009al), or even evolved stars ^^at planet formation is not rare. 
(Johnsonl l2007l: |Lovis & Mayor 2007). The sample of detected 

and characterized planets thus becomes large enough to per- Observational evidences regai'ding close-in planets lead to 

form robust statistical analysis of the population and test plan- favor a formation by the co r e-accreti on mechanism (hereafter 

etary formation theories. In that sense, the planet occuiTence CA, e.g. |Pollack et alj|l996t). |Sousa et alj (|20il) show that the 

frequency has been determined for giant and telluric planets, presence of close-in giant planets is correlated with the metal- 

licity of their host stars. This correlation is al so related to, if 

* Based on observations collected at the European Organization for Pl^nets orbit within 3 AU, their host-star mass Lovis & Mayor] 

Astronomical Research in the Southern Hemisphere, Chile, ESO : (2007); Johnson et al. (2010); Bowler et al. (2010). Another cor- 

runs 084.C-0396A, 085.C-0675A, 085.C-0277B, 087.C-0292A, 087.C- relation but between the content in heavy elements of the planets 

0450B, 088.C-0085A, 089.C-0149A. and the metallicity of their parent star has also been found (e.g. 
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Table 1. Frequency of giant planets reported by various surveys around the full spectral type range. 



Sample 


ieciinique 


Sep. range 


Mass range 


Frequency 


Distribution 


Reference 






^ A T T\ 

(AU) 


(Mj) 


(%) 


(AU) 




102 M 


RV 


< 1 


< 3 


1-5 


observed 


Bonfils etal. ('2013') 


822 FGK 


RV 


<5 


>0.3 


14 


observed 


Mavoretal. (2011) 


31 old- A 


RV 


0.1 - 3 


0.5 - 14 


9-14 


observed 


Johnson etal. (2010) 


585 F-M 


RV 


< 3 


0.3 - 10 


10.5 


observed 


Cummins et al. (2008) 


42 AF 


AO 


5-320 


3-14 


5.9-18.8 


flat / Cu08° 


Visanetal. (2012) 


85 F-M 


AO 


5-500 


< 100 


< 10 


flat + Gl'' 


^anson et al. (2012) 


15 BA 


AO 


< 300 


< 100 


< 32 


flat + Gl 


Janson et al. (201 1) 


118 F-M 


AO 


25 - 856 


> 4 


< 20 


flat / Cu08 


Nielsen & Close (2010)" 


88 B-M 


AO 


>40 


5-13 


< 10 


Cu08 


Chauvin et al. (2010) 


88 FGKM 


AO 


50 - 250 


0.5-13 


< 9.3 


power laws in m and a'' 


Lafreniere et al. (2007) 


22 GKM 


AO 


> 30 


> 3 


< 5 


power laws in in and a'' 


Kasoer et al. (2007) 



Notes. Results on the frequency of giant planets are reported at 68 % confidence level, except the ones from 'Janson et al.l ( |2011| . |20T21) which are 
stated at 99 % confidence level. Cummina et al. (2008) They infer the planet population from boundaries in a planet m ass-semi ma j or axi s 
grid considering disk instability model. They performed their analysis using results from surveys of Masciadr i et all ( l200^ : lBiller et al.l ( 120071) : 
Lafreniere et al. (2007). They infer the planet population from power laws distributions with different coefficients to m and a as the ones in 
Cumming et al.l ( l2008h . Please be referred to the publication for details. 



iGuillot et al.ll2006HMiller & Fortnevll201 ih . All are insights that 
CA is operating at short orbits. According to this scenario, the 
first steps of the growth of giant gaseous planets are identical to 
those of rocky planets. The dust settling towards the mid plane 
of the protoplanetary disk leads to formation of larger and larger 
aggregates through coagulation up to meter-sized planetesimals. 
These cores grow up then through collisions with other bod - 
ies until they reach a critical mass of 10 Me dM 

izun 

Their gravitational potentials being high enough, they trigger 
runaway gas accretion and become giant planets. However, such 
scenario requires high surface density of solids into the disk 
to provide enough material to form the planet core and a large 
amount of gas. Large ga seous planets are not expected t o form 
in situ below the ice Un e. lLin et all (ll996h : lAiibert et al.l (|2004 : 
iMordasini et al.l (l2009h refine the model with inward migration 
to explain the large amount of giant planets orbiting very close 
to their parent stars. 

At wider (> 30 AU) orbits, the situation is very dif- 
ferent since this core accretion mechanism has difficu lties to 
form giant planets (Bolevetal. 2009; Dodson-Robinso n et al.l 
120091) . The timescales to form massive cores become longer 
than the gas dispersal ones and the disk surface density too 
low. Additional outward migration mec hanisms must be invok ed 
(corotation torque in radiative disks, Ki ev & NelsonI 1201 2l or 
planet-planet scattering, Cridaet al. 2009). Alternatively, cloud 
fragme ntation can form objec ts down to the planetary mass 
regime ( iWhitworth et al.ll2007l) and is a solid alternative to ex- 
plain the existence of very wide orbits substellar companions. 
Finally, disk fragmentation also called gravitational instability 
(Gl, Cameron 1978; Stamatellos & Whitworth 2009) remains an 
attractive mechanism for the formation of massive giant planets 
beyond 10 to 20 AU. According to this scenario, a protoplane- 
tary disk becomes unstable if cool enough leading to the exci- 
tation of global instability modes, i.e. spiral arms. Due to their 
self-gravity, these arms can break up into clumps of gas and dust 
which are the precursor for giant planets. 

Understanding how efficient are these different mechanisms 
as a function of the stellar mass, the semi-major axis, and the 
disk properties are the key points to fully understand the for- 
mation of giant planets. Understanding how giant planets form 



and interacts with their environment is crucial as they will ulti- 
mately shape the planetary system's architecture, drive the tel- 
luric planet's formation, and the possible existence of conditions 
favorable to Life. 

The presence of massive dusty disks around young stars, 
like HR 8799 and /3 Pictoris, might be a good indicator of the 
presence of exoplanetary systems recently formed (Rhee et aT] 
1200 7). Observations at several wavelengths revealed asymme- 
try structures, linglike sometimes, or even warps which could 
arise from gravitat ional perturbat i ons imposed by one or more 
gia nt planet (e.g . iMouillet et"d] Il997t lAugereau et al.l 120011: 
Kalas et al. l l2005h. Thanks to improvement of direct imaging 
(DI) technique with ground-based adaptive optics systems (AO) 
or space telescopes, a few planetary mass objects and low 
mass brown dwarfs have been detected since the first one by 
Cha uvin et al. 1(12004'). One also has to refer to the breakthrough 
discoveries of giant planets between 8 an d 68 AU around 
young, nearby, and d u sty early-type stars (iKalas et al.1 20081: 
Lagrange et akl l2009bl: iMarois et al.l l2008l 120101: ICarson et al.l 



20121) . Direct imaging is the only viable technique to probe for 



planets at large separations but detecting planets need to over- 
come the difficulties due to the angular proximity and the high 
contrast involved. 

Nevertheless, numerous large direct imaging surveys to 
detect giant pla net companions have repotted null detection 



(Masciadri et aLll2005t iBiller et al.ll2007l: lUa^ 
Ehrenreich et al.^'2 ()10l: IChauvin et al.ll2010r 



reniere et al. 



Janson et al.l 



2007 



2011 



Delorme et al. 2012|). nevertheless this allowed to set upper lim- 
its to the occurrence of giant planets. Table [T] reports the sta- 
tistical results of several direct imaging surveys, as a function 
of the sample, separation and mass ranges , and planet distri- 
bution. All surveys previous to the one of'Vigan et all (1201 2l) 
derive upper limits to the occurrence of giant planets, usually 
more massive than 1 - 3Mj between few to hundreds of AUs. 
They find that less than 10 - 20 % of any star harbor at least 
one giant planet if the distribution is flat or similar to the RV 
one, taking into account all the assumptions beyond this results. 
Janson et al. (201 1, 2012) include in the planet distribution lim- 
itations if giant planets form via Gl. They show that the occur- 
rence of planets might be higher for high mass stars than for 
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solar- type s tars but GI is still a rare formation channel. On the 
other hand.'Vigan et alJ (1201 2*) takes into account two planetary 
system detections among a volume-limited set of 42 A-type stars 
to derive lower limits for the first time. It comes out that the fre- 
quency of jovian and massive giant planets is higher than 5.9 % 
around A-F stars. However, all these surveys suggest a decreas- 
ing distribution of planets with increasing separations, which 
counterbalances the RV trend. 

In this paper we report the results of a deep direct imaging 
survey of 59 young, nearby, and dusty stars aimed at detecting 
giant planets on wide orbits performed between 2009 and 2012. 
The selection of the target sample and the observations are de- 
tailed in Section |2] In Section [3] we describe the data reduction 
and analysis to derive the relative astrometry and photometry of 
companion candidates, and the detection limits in terms of con- 
trast. Section|4]is then dedicated to the main results of the survey, 
including the discovery of new visual binaries, the characteriza- 
tion of known substellar companions, and the detection perfor- 
mances. Finally, we present in Section |5] the statistical analysis 
over two special samples: A-F type stars and A-F dusty stars, 
from which we constrain the frequency of planets based on dif- 
ferent formation mechanisms or planet population hypotheses. 

2. Target sample and observations 

2.1. Target selection 

The target stars were selected to optimize our chance of planet 
detection according to : 

- Distance : with a given angular resolution limited by the 
telescope's diffraction limit, the star's proximity enables to 
access closer physical separations and fainter giant planets 
when background-limited. We therefore limited the volume 
of our sample to stars closer than 150 pc, even closer than 
100 pc for 94 % of them (see Figure [1] top left panel). 

- Observability and magnitude : stars were selected according 
to : 1) their declination {S < 25 deg) for being observable 
from the southern hemisphere, 2) their K-band brightness (K 
< 7 mag) to ensure optimal AO corrections, 3) for being 
single to avoid degradation of the AO performances and 4) 
for being never observed in deep ima ging (see section [2}. 

- Age : evolutionary models ( Baraffe et a l. 2003: Marl ev et al] 
12007) predict that giant planets are intrinsically more lumi- 
nous at young ages and become fainter with time. Therefore, 
for a given detection threshold, observing younger stars is 
sensitive to lower mass planets. Our sample selection is 
based on recent publications on associations (AB Dor, /3 Pic, 
Her/Lyr, Argus, Tuc/Ho r, Columba, Upper C en/L upus) from 
IZuckerm an et al.1 (1201 ih . iTorres et al.l (l2008b and iRhee et all 
( 12007 ). Indeed, stars belonging to these moving groups share 
common kinematic properties and ages. These parameters 
are measured from spectroscopy, astrometry, and photome- 
try (optical and X-rays). 64 % of the selected stars belongs 
to young and nearby moving groups. 88% of the targets are 
younger than 200 Myr and even 62 % younger than 70 Myr 
(see Figure[Tl top right panel). 

- Spectral Type : recent imaged giant planets have been de- 
tected around the intermediate-mass HR8799, Fomalhaut, 
and /3 Pictoris stars with separations from 8 to 110 AU. More 
massive stars imply more massive disks, which potentially 
allow the formation of more massive planets. We therefore 
have biased 79 % of our sample towards spectral types A 
and F (see Figure[T] bottom panel). 



Table 4. Sample definitions 



Name 


Number 


SpT 


Representative rate 

(%) 


Survey 


59 


B-M 




A-F sample 


37 


A-F 


56 


A-F dusty sample 


29 


A-F 


72 



- Dust : dusty debris disks around pre- and main-sequence 
stars might be signposts for the existence of planetesimals 
and exoplanets (see a review in Krivov 2010). 76 % of our 
sample are star with large infrared excess at 24 and/or 70 fim 
(IRAS, ISO and Spitzer/MlPS), indicative of the emission of 
cold dust. The remaining stars have no reported excess in the 
literature. 

The name, coordinates, galactic latitude (b), proper motions 
(ju), spectral type (SpT), distance (d), K magnitude, and age of 
the target stars are listed in Table |2] together with the reference 
for the age determination and the moving group they belong to if 
they dqj Figure [T] summarizes the main properties of the target 
stars. Briefly, the sample consists of 59 B- to M-type stars whose 
the median star would be a F-type at distance of 40 pc with an 
age of 30 Myr, a K-magnitude of 5.5, and an apparent proper 
motion of 85 mas/yr 

To analyze an homogeneous and volume-limited sample, we 
perform the statistical study on stars which have 1) d < 65 pc, 
2) age < 100 Myr 3) dec < 25 deg, and 4) Spectral type - A 
or F. We get set of 68 young, nearby, and southern A-F stars 
from the literature. 33 stars in our survey fulfill these criteria 
(flagged with a ★ symbol in Table|2]l, i.e. 48 % of completeness. 



from a previous survey (IVigan et 



ing observ( 

in i2or2iP i 



(see Table O, thus 
reaching a representative rate of 56 %. We will then refer to this 
sample of stars as the A-F sample. 

Moreover, we also aim at constraining the formation mech- 
anism and rate of giant planets around /3 Pictoris analogs. A 
sub-set of stars is extracted from the A-F homogeneous sample 
by considering an IR excess at 24 and/or 70 urn (from the 
same references as for our survey plus M izusawa et al. 2012[ 
'Rebull et al.' '2008'; 'Hillenbrand etaP 120081: iSu et alJ llOm . 
Among the full se t of 68 stars , 39 are dusty. Our survey plus 
1 stars from Vie a n et al.l (l2012h . i.e. 28 stars, reach a complete 
level of 72 %. We will then refer to this sample as the A-F dusty 
sample. 

Table m summarizes all different samples. 



2.2. Observing strategy 

The survey was conducted betwee n 2009 and 2012 w ith the 
NAOS adaptive optics instrument (iRousset et al.l l2003b com- 



' We attempt to derive the age in a homogeneous way. If the star 
belongs to a moving group, the age of that group is adopted for this star. 
If the star does not belong to a known association, then we ensured that 
the age determination was done on similar way than for the membership 
identification, i.e. the galactic space motions UVW, the Li /16708 A line 
equivalent width or the X-ray emission. 

- These additional stars have been also observed using ADI tech- 
niques but with the Ks filter on VLT/NaCo. 
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Table 2. Sample of young, nearby, and dusty stars observed during our VLT/NaCo thermal and angular differential imaging survey. 



Name 


a 6 


b 


Ha C0S((5) ytis 


d 


SpT 


K 


excess ? age Ref. 


HIP HD 


(J2000) (J2000) 


(deg) 


(mas.yr"') (mas.yr"') 


(PC) 




(mag) 


(Myr) 


AB Doradus 



6276 




01 20 32 


-1128 03 


-72.9 


110.69 


-138.85 


35.06 


G9V 


6.55 


y 


70 


★18859 


25457 


04 02 37 


-00 16 08 


-36.9 


149.04 


-253.02 


18.83 


F7V 


4.18 


y 


70 


30314 


45270 


06 22 31 


-60 13 08 


-26.8 


-11.22 


64.17 


23.49 


GIV 


5.04 


y 


70 


★93580 


177178 


19 03 32 


01 49 08 


-1.86 


23.71 


-68.65 


55.19 


A4V 


5.32 


n 


70 


95347 


181869 


19 23 53 


-40 36 56 


-23.09 


32.67 


-120.81 


52.08 


B8V 


4.20 


n 


70 


109268 


209952 


22 08 13 


-46 57 38 


-52.47 


127.6 


-147.91 


31.09 


B6V 


2.02 


n 


70 


★ 115738 


220825 


23 26 55 


01 15 21 


-55.08 


85.6 


-94.43 


49.7 


AO 


4.90 


y 


70 


★117452 


223352 


23 48 55 


-28 07 48 


-76.13 


100.03 


-104.04 


43.99 


AOV 


4.53 


y 


70 



(3 Picloris 



★ 11360 


15115 


02 26 16 


+06 17 34 


-49.5 


86.09 


★21547 


29391 


04 37 36 


-02 28 24 


-30.7 


43.32 


★25486 


35850 


05 27 05 


-11 54 03 


-24.0 


17.55 


★27321 


39060 


05 47 17 


-51 03 59 


-30.6 


4.65 


★27288 


38678 


05 46 57 


-14 49 19 


-20.8 


14.84 


★79881 


146624 


16 18 18 


-28 36 50 


-1-15.4 


-33.79 


★88399 


164249 


18 03 03 


-51 38 03 


-14.0 


74.02 


★92024 


172555 


18 45 27 


-64 52 15 


-23.8 


32.67 


★95261 


181296 


19 22 51 


-54 25 26 


-26.2 


25.57 


★95270 


181327 


19 22 59 


-54 32 16 


-26.2 


23.84 


102409 


197481 


20 45 09 


-31 20 24 


-36.8 


280.37 



-50.13 


44.78 


F2 


5.86 


y 


12 


10 


-64.23 


29.76 


FOV 


4.54 


n 


12 


2 


-50.23 


27.04 


F7V 


4.93 


y 


12 


2 


83.1 


19.4 


A6V 


3.53 


y 


12 


2 


-1.18 


21.52 


A21V/V 


3.29 


y 


12 


13 


-100.59 


43.05 


AOV 


4.74 


n 


12 


2 


-86.46 


48.14 


F5V 


5.91 


y 


12 


2 


-148.72 


29.23 


A7V 


4.30 


y 


12 


2 


-82.71 


48.22 


AOVn 


5.01 


y 


12 


2 


-81.77 


50.58 


F6V 


5.91 


y 


12 


2 


-360.09 


9.94 


MIV 


4.53 


y 


12 


2 



Tucana-Horologium / Columba 



★1134 


984 


00 14 10 


-07 11 56 


-66.36 


102.84 


-66.51 


46.17 


F7V 


6.07 


n 


30 


★2578 


3003 


00 32 44 


-63 01 53 


-53.9 


86.15 


-49.85 


46.47 


AOV 


4.99 


y 


30 


7805 


10472 


01 40 24 


-60 59 57 


-55.1 


61.94 


-10.56 


67.25 


F21VA^ 


6.63 


y 


30 


★9685 


12894 


02 04 35 


-54 52 54 


-59.2 


75.74 


-25.05 


47.76 


F4V 


5.45 


n 


30 


10602 


14228 


02 16 31 


-51 30 44 


-22.2 


90.75 


-21.9 


47.48 


BOV 


4.13 


n 


30 


12394 


16978 


02 39 35 


-68 16 01 


-45.8 


87.4 


0.56 


47.01 


B9 


4.25 


n 


30 


16449 


21997 


03 31 54 


-25 36 51 


-54.1 


53.46 


-14.98 


73.80 


A3IVA' 


6.10 


y 


30 


★22295 


32195 


04 48 05 


-80 46 45 


-31.5 


46.66 


41.3 


61.01 


F7V 


6.87 


y 


30 


★26453 


37484 


05 37 40 


-28 37 35 


-27.8 


24.29 


-4.06 


56.79 


F3V 


6.28 


y 


30 


26966 


38206 


05 43 22 


-18 33 27 


-23.1 


18.45 


-13.2 6 


69.20 


AOV 


6.92 


y 


30 


★30030 


43989 


06 19 08 


-03 26 20 


-8.8 


10.65 


-42.47 


49.75 


F9V 


6.55 


y 


30 


30034 


44627 


06 19 13 


-58 03 16 


-26.9 


14.13 


45.21 


45.52 


KIV 


6.98 


y 


30 


★107947 


207575 


21 52 10 


-62 03 08 


-44.3 


43.57 


-91.84 


45.09 


F6V 


6.03 


y 


30 


★114189 


218396 


23 07 29 


+21 08 03 


-35.6 


107.93 


-49.63 


39.40 


FOV 


5.24 


y 


30 


★118121 


"'24392 


23 57 35 


-64 17 53 


-51.8 


78.86 


-61.1 4 


8.71 


AlV 


4.82 


n 


30 



Argus 



★67945 08 09 39 -20 13 50 +7.0 -38.6 25.8 63.98 FOV 7.15 n 40 3 

★57632 102647 11 49 04 +14 34 19 +70.8 -497.68 -114.67 11.00 A3V 1.88 y 40 3 



Hercules-Lyra 



544 166 00 06 37 +29 01 19 -32.8 379.94 -178.34 13.70 KOV 4.31 y 200 4 

7576 10008 01 37 35 -06 45 37 -66.9 170.99 -97.73 23.61 G5V 5.70 y 200 4 



Upper Centaurus-Lupus 



78092 142527 15 56 42 -42 19 01 -11.19 -24.46 145. F6nie 4.98 y 5 6 



Other 



682 


377 


00 08 26 


+06 37 00 


+20.6 


88.02 


★7345 


9692 


01 34 38 


-15 40 35 


-74.8 


94.84 


7978 


10647 


01 42 29 


-53 44 26 


-61.7 


166.97 


★ 13141 


17848 


02 49 01 


-62 48 24 


-49.5 


94.53 


18437 


24966 


03 56 29 


-38 57 44 


-49.9 


29.46 


22226 


30447 


04 46 50 


-26 18 09 


-37.9 


34.34 


★22845 


31295 


04 54 54 


+10 09 03 


-20.3 


41.49 


34276 


54341 


07 06 21 


-43 36 39 


-15.8 


5.8 


38160 


64185 


07 49 13 


-60 17 03 


-16.6 


-37.41 


★41307 


71155 


08 25 40 


-03 54 23 


+ 18.9 


-66.43 



-1.31 


39.08 


G2V 


6.12 


y 


30 


7 


-3.14 


59.4 


AlV 


5.46 


y 


20 


5 


-106.71 


17.35 


F9V 


4.30 


y 


300 


5 


29.02 


50.68 


A2V 


5.97 


y 


100 


5 


0.1 


105.82 


AOV 


6.86 


y 


10 


5 


-4.63 


78.125 


F3V 


6.89 


y 


100 


5 


-128.73 


35.66 


AOV 


4.41 


y 


100 


5 


13.2 


102.35 


AOV 


6.48 


y 


10 


5 


140.08 


34.94 


F4V 


4.74 


n 


200 


8 


-23.41 


37.51 


AlV 


4.08 


y 


100 


5 
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Table 2. Continued 



Name 


a 


6 


b 


fia COS(6) 




d 


SpT 


K 


excess ? 


age 


Ref. 


HIP 


HD 


(J2000) 


(J2000) 


(deg) 


(mas.yr"') 


(mas.yr"') 


(pc) 




(mag) 




(Myr) 




JOJZ4 




10 57 03 


-68 40 02 


-8.1 


-41.41 


12.47 


90.42 


A8III 


6.79 


y 


50 


5 


59315 


105690 


12 10 07 


-49 10 50 


+ 13.1 


-149.21 


-61.81 


37.84 


G5V 


6.05 


n 


100 


9 


76736 


138965 


15 40 12 


-70 13 40 


-11.9 


-40.63 


-55.31 


78.49 


A3V 


6.27 


y 


20 


5 


★86305 


159492 


17 38 06 


-54 30 02 


-12.0 


-51.04 


-149.89 


44.56 


A7IV 


4.78 


y 


50 


11 


★99273 


191089 


20 09 05 


-26 13 27 


-27.8 


39.17 


-68.25 


52.22 


F5V 


6.08 


y 


30 


5 


★101800 


196544 


20 37 49 


+ 11 22 40 


-17.5 


39.15 


-8.26 


7.94 


AlIV 


5.30 


y 


30 


5 


108809 


209253 


22 02 33 


-32 08 00 


-53.2 


-19.41 


23.88 


30.13 


F6.5V 


5.38 


y 


200 


5 


1 14046 


217987 


23 05 47 


-35 51 23 


-66.0 


6767.26 


1326.66 


3.29 


M2V 


3.46 


n 


100 


12 




219498 


23 16 05 


+22 10 02 


-35.6 


79.7 


-29.4 


150.0 


G5 


7.38 


y 


300 


7 


★116431 


221853 


23 35 36 


+08 22 57 


-50.0 


65.37 


-40.79 


68.45 


FO 


6.40 


y 


100 


5 



Notes. Stars with the ★ sym bol are used for the statistical analysis. Star parameters (a, S, b, yu„ cos((5), yu^ and d) are extracted from the Hipparcos 
catalog ('van Leeuwen'2007). Unit of right ascension are hours, minutes, and seconds ; units of declination are degrees, arcminutes, and arcseconds. 
For HD 219498, the distance is extracted from Roccatagliata et al. (2009). The K magnitudes ar e extra c ted from the 2MASS catalog (Cutri et ah] 
12003 V The IR excess at 24 and/or 70 jim are extracted from Zuckerman et al. (201 1); Kains et al.l (l201 ll) : iMorales et alJ l l2011h : lRhee et alj ( 2007). 
The age of HIP 93580 is still debated in Zuckerman et a l. (201 1) due to discordant kinematics. The binarity (if physical) may have some impact 
on both proper motions and RVs and thus on the membership to AB Dor. 
The age references are the following : 



Refer ences. ( 1 ) IZuckerman et al .l (!20 1 ll) : (2) IZuckerman et al.l ( |200Tb ; (3) [Torres et"al1 (120081) ; (4)lL6pez-Santiag o et al.l ( 120061) ; (5)lRhee et al 
'20071) ; (6) See discusion in|Rameau et al. ( 2012) ; (7) Hillenbrand et al. (2008) ; (8) Zuckerman et al. (2006) ; (9) See for instancelChauvin et al. 



20ld) ; (lO) ISchliederet al .1120121) ; (ll) ,Song et al..(i2001i) ; (12) See discusion in iDelorme et al., ( i2012i) ; (13) ,Nakajima & Morind ( 120121) 



Table 3. Additional young, nearby and southern A-F stars for the statistical analysis 



HIP 



Name 
HD 



a 

(J2000) 



6 

(J2000) 



b 

(deg) 



fi„ cos((5) 
(mas.yr"') 



(mas.yr"') 



d 

(pc) 



SpT 



K 

(mag) 



excess ? 



age 
(Myr) 



Ref. 



★ 12413 


16754 


02 39 48 


-42 53 30 


-63.0 


88.20 


-17.82 


39.8 


AlV 


4.46 


n 


30 


1 


★ 14551 


19545 


03 07 51 


-27 49 52 


-59.8 


66.26 


-19.09 


54.6 


A5V 


5.77 


n 


30 


1 


★26309 


37286 


05 36 10 


-28 42 29 


-28.1 


25.80 


-3.04 


56.6 


A2III 


5.86 


y 


30 


1 


★61468 


109536 


12 35 45 


-41 01 19 


21.8 


-107.09 


0.63 


35.5 


A7V 


4.57 


n 


100 


1 



Notes. The reference corresponds to the publication the star is extracted from : (l) IVigan et al.l(l2012h . 



bined to the CONICA near-infrared camera (iLenzen et al 
NaCo is mounted at a Nasmyth focus of one of the 8.2 m ESO 
Very Large Telescopes. It provides diffraction-limited images 
on a 1024 x 1024 pixel Aladdin 3 InSb array. Data were ac- 
quired using the L27 camera, which provides a spatial sam- 
pling of ^ 27.1 mas/pixel and a field of view (FoV hereafter) 
of 28 " X 28 ". In order to maximize our chance of detection, 
we used thermal-infrared imaging with the broadband L ' fil- 
ter (Ac - 3.8 fim, AA - 0.62 fim) since it is optimal to detect 
young and warm massive planets with a peak of emission around 
3-4 jUOT. 

NaCo was used in pupil-tracking mode to reduce instrumen- 
tal speckles that limit the detection performances at inner angles, 
typically between 0. 1 " and 2.0". This mode provides rotation of 
the F oV to use angular differential imaging (ADI, iMarois et all 
I2OO 6). The pupil stabilization is a key element for the second 
generation instruments GPI and VLT/SPHERE. Note that NaCo 
suffered from a drift of the star with time (few pixels/hours de- 
pending on the ele vation) associated t o the pupil-tracking mode 
until October 201 1 dGirard et al.ll2012l) . Higher performance was 
obtained after the correction of the drift. To optimize the image 



selection and data post-processing, we recorded short individual 
exposures coupled to the windowing mode of 512 x 514 pixels 
(reduced FoV of ^ 14" x 14"). The use of the dithering pat- 
tern combined to the cube mode also ensure accurate sky and 
instrumental background removal. The detector integration time 
(DIT) was set to 0.2 s to limit the background contribution to the 
science images. 

Each observing sequence lasted around 90 minutes, includ- 
ing telescope pointing and overheads. It started with a sequence 
of short unsaturated exposures at five dither positions with the 
neutral density filter NDio„g (transmission of 1.17 %). This al- 
lowed the estimation of the stellar point spread function (PSF) 
and served as photometric calibrator. Then, saturated science 
images were acquired with a four dithering pattern every two 
DITxNDIT exposures with NDIT= 100 stored into a datacube 
and this was repeated over more than 100 times to provide suf- 
ficient FoV rotation for a given star. The PSF core was saturat- 
ing the detector over a ^ 5 pixel-wide area. Twilight flat-fields 
were also acquired. For some target stars, second epoch data on 
NaCo were acquired with the same observing strategy. Finally, 
9i Ori C field was observed as astrometric calibrator for each 
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d(pc) 



Age (Myr) 




FO GO KO MO 
Spectral type 

Fig. 1. Graphical summaries of the main properties on the target stars. Top raw : Left : Histogram of distances with 10 pc bin. 
Right : Histogram of ages of known members of considered moving groups (AB Dor - 70 Myr, y6 Pictoris - 12 Myr, Tuc/Hor & Col 
- 30 Myr, Her/Lyr - 200 Myr, Argus - 40 Myr, Upper Cen/Lup -5/10 Myr) and additional stars with 1 Myr bin. Central raw : Left 
: Histogram of K-band magnitude with 0.5 mag bin. Right : Histogram of apparent proper motions with 25 mas/yr bin. Bottom 
raw : Histogram of spectral types. 



observin g run. The same set of stars ori ginally observed with 
HST bv lMcCaughrean & StauffeJ (11994 (TYC058, 057, 054, 
034 an 026) were imaged with the same set-up (L' with the L27 
camera). The mean platescale and tiTie North orientation were 
measured and reported on Table |5] 

2.3. Observing conditions 

Observations in period 84 and 85 were done in visitor mode as 
the pupil stabilization mode was not offered in service mode. 
Observations in period 87, 88 and 89 were then completed in 
service mode to benefit from optimal atmospheric conditions. 
A summary of the observing conditions is reported on Figure 



|2] showing histograms of explored parallactic angle ranges, air- 
mass, as well as the atmospheric conditions : seeing and coher- 
ence time To. Note that NAOS corrects the atmospheric turbu- 
lences for bright stars when tq remains longer than 2 ms (63 % 
of the observations). When tq decreases, the image quality and 
precision for astrometric and photometric measurements are de- 
graded. The observations were however conducted under good 
conditions since the median seeing is 0.9 ", the median tq is 
3.2 ms and the median airmass is 1.15. Finally, 72% of the stars 
were observed with a parallactic angle exploration larger than 
20 deg. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 2 4 6 8 10 

Seeing (arcsec) TauO (ms) 




Airmass Explored parang, (deg) 

Fig. 2. Graphical summaries of the log of the observations of the target young, nearby dusty stars observed with VLT/NaCo between 
2009 and 2012. Top left : Histogram of median seeing (image quality seen by the active optics sensor) with 0.2" bin. The DIMM 
seeing has been over plotted with dashed columns. Top left : Histogram of median tq with 0.5 ms bin as calculated by NAOS. 
Bottom left : Histogram of median airmass with 0.1 bin. Bottom right : Histogram of explored parallactic angles with 10 deg bin. 



Table 5. Mean platescale and true north orientation measured on 
6*1 Ori C field with NaCo and L'/L27 set-up. 



ESO program UT-date 



Platescale 
(mas) 



True north 
(deg) 



084.C-0396A 


11/24/2009 


27.09 ± 0.02 


-0.08 


±0.01 


085.C-0675A 


07/27/2010 


27.12 ±0.03 


-0.36 


±0.05 


085.C-0277B 


09/28/2010 


27.11 ±0.04 


-0.36 


±0.11 


087.C-0292A 


12/18/2011 


27.10 ±0.03 


-0.60 


±0.01 


087.C-0450B 


12/08/2011 


27.16 ±0.08 


-0.52 


±0.07 


088.C-0885A 


02/19/2011 


27.10 ±0.03 


-0.38 


±0.03 


089.C-0149A 


08/24/2012 


27.11 ±0.02 


-0.41 


±0.07 



3. Data reduction and analysis 

3.1. Unsaturated images 

The unsaturated dithered exposures of each star were p roces sed 
with the Eclipse software developed bv lDevillaJ (1 1997b : bad pix- 
els removal, sky subtraction constructed as the median of the im- 
ages followed by flat-fielding were applied to data; the final PSF 
image was then obtained by shifting and median combining the 
images. 



3.2. Saturated angular differential images 

The reduction of the ADI saturated dithered datacubes was per- 
formed with the dedicated pipeline developed at the Institut 
de Planetologie et d'Astrophysique de Grenoble (IPAG). This 
p ipeline has been inte n sively used and gave p robing results 
: 'Lagrange et al.' ('2010'); 'B onnefov etal] (1201 ll) ; IChauvin et al] 
(2012); Delorme et al. (2012i); Lagrange et al. (2012); Rameau 
et al. (2012). We describe the main steps in the following. 

Getting twilight flats allowed us to achieve optimal flat- 
fielding and bad pixel identification. We used the Eclipse soft- 
ware to extract those calibrations frames. The raw data were then 
divided by the flat-field and removed for the bad and hot pixels 
through interpolation of the closest neighbor pixels. Sky estima- 
tion was performed by taking the median of the 400 s closest 
in time dithered exposures within a cube and then subtracted 
to each frame. Frames with low quality were removed from the 
cubes following a selection based on cube statistics such as the 
flux maximum, the total flux, and the encircled energy in each 
frame in an annulus outside the saturated pixels. Poor quality 
frames due to degraded atmospheric conditions were rejected 
(typically less than 10% of the complete observing sequence, 
see Girard et al. 2010 for the cube advantages). The good-quality 
frames were recentered to a common central position us ing the 
Eclipse software for the shift and Moffat profile fitting (IMoffatl 
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Il969h on the PSF wings for the registration of the central star 
position. We ended up with good-quality cleaned and recentered 
images within a single master cube associated with their paral- 
lactic angle values. A visual inspection was done to check the 
quality of the final frames. 

Subsequent steps were the estimation and subtraction of the 
stellar halo for each image then derotation and stacking of the 
residuals. The most critical one is the estimation of the stellar 
halo which drives the level of the residuals. We applied different 
ADI algorithms to optimize the detection performances and to 
identify associated biases. Since the quasi-static speckles limit 
the performances on the inner part of the Fo V, we performed the 
ADI reduction onto reduced frames, typically 200 x 200 pixels. 
We recall here the difference between the four ADI procedures : 

- in classic ADI (cADI, 'Marois et al."2006'), the stellar halo 
is estimated as the median of all individual reduced images 
and then subtracted to each frame. The residuals are then 
median-combined after the derot ation; 

- in smart ADI (sADI, Lagran ge et all |2010|) . the PSF- 
reference for one image is estimated as the median of the 
n closest-in-time frames for which the FoV has rotated more 
than QfX FWHM at a given separation. Each PSF-reference 
is then subtracted to each frame and the residuals are mean- 
stacked after the derotation; We chose a PSF-depth of n = 10 
frames for the PSF-building satisfying a separation criteria of 
a = 1. FWHM at a ra dius of 1.3"; 

- the radial ADI (rADI. lMarois et al.l2006t) procedure is an ex- 
tension of the sADI where the n frames with a given rotation 
used for the stellar halo building are selected according to 
each separation. The PSF-depth and the a coefficient were 
chosen as for sADI (n - 10 and a = 1. FWHM). The radial 
extent of the PSF-building zone is Ar = 1 .4 FWHM below 
1 .6 " and 3 FWHM beyon± 

- in the LOCI approach dLafreniere et all l2007b , the PSF- 
reference is estimated for each frame and each location 
within this frame. Linear combinations of all data are com- 
puted so as to minimize the residuals into an optimization 
zone, which is much bigger than the subtraction zone to 
avoid the self-removal of point-like sources. We considered 
here a radial extent of the subtraction zone Ar = 0.9 FWHM 
below 1 .6 " and 3 beyond; a radial to azimuthal width ra- 
tio was set to ^ = 1 ; a standard surface of the optimization 
zone was A^^ - 300 PSF cores; the separation criteria of 
A^<5 = 1 FWHM. 

All the target stars were processed in a homogeneous way 
using similar set of parameters. It appears that when the PSF re- 
mained very stable during a sequence (i.e. tq > 4 ms), advanced 
ADI techniques do not strongly enhance the performance. 

ADI algorithms are not the best performant tools for 
background-limited regions as the PSF-subtraction process add 
noise. We thus processed the data within the full window (i.e. 
512 X 514 pixels with the dithering pattern) with what we called 
the non-ADI (nADI) procedure. It consists in 1) computing an 
azimuthal average of each frame within 1 -pixel wide annulus, 2) 
circularizing the estimated radial profile 3) subtracted the given 
profile to each frame and then 4) derotating and mean-stacking 
the residuals. nADI by-products can help to distinguish some 
ADI artifacts from real features as well. 

For each star, a visual inspection of the five residual maps 
was done to look for candidate companions (CC). 



3.3. Relative photometry and astrometry 

Depending on the separation and the flux of the detected CC, 
different techniques were used to retrieve the relative photometry 
and astrometry with their uncertainties: 

- for brig ht visu al binaries, we used the deconvolution algo- 
rithm of IVeran & Rigaut ( 1998); 

- for CCs detected in background-limited regions (in nADI fi- 
nal images), the relative photometry and astrometry were ob- 
tained usi ng a 2D moffat fittin g and classical aperture pho- 
tometry ( Chauvin et al.ll20Toh . The main limitation of this 
technique remains the background subtraction which affects 
the level of residuals; 

- for speckle limit obj ects, fake planet s were i njected follow- 
ing th e approach of iBonnefov et al.l (l201lh : Ichauvin et al.l 
( I2OI2 !) with the scaled PSF-reference at the separation of the 
CC but at different position angles. The injections were done 
into the cleaned mastercubes which were processed with the 
same setup. We then measured the position and the flux of 
the fake planets which minimized the difference with the real 
CC. The related uncertainties associated with this method 
were also estimated using the various set of fake planets in- 
jected at different position angles. 

In both algorithms, the main error for the relative astrom- 
etry is the actual center position of the saturated PSF (up to 
0.5 pixel). Other sources of errors come from the Moffat fit- 
ting, the self-subtraction, the residual noise, or the PSF shape. 
The reader can refer for more details to the dedicated analy- 
sis o n uncertainties on C C astrometry using VLT/NaCo ADI 
data (IChauvin et al ] |20T1 . For a CC observed at several epochs 
(follow-up or archive), we investigated its status (background 
source or comoving object) by determining its probability to be 
a stationary background object, assuming no orbital motion. This 
approach is the same as in Chauvin et al. (2005) by comparing 
the relative positions in a and 6 from the parent-star between 
the two epochs, from the expected evolution of positions of a 
background object, given the proper and parallactic motions and 
associated error bars. 

3.4. Detection limits 

The detection performances reached by our survey were es- 
timated by computing 2D detection limit maps, for each tar- 
get star, at 5cr in terms of L ' contrast with respect to the pri- 
mary. For each set of residual maps for each target, we com- 
puted the pixel-to-pixel noise within a sliding box of 1 .5 x 1 .5 
FWHM. The second step was to estimate the flux loss due to self- 
subtraction by the ADI processing. We created free-noise-cubes 
with bright fake planets (100 ADU) at three positions, i.e. 0, 120 
and 240 deg, each 20 pixels from the star, with the same FoV 
rotation as real datacubes for each star then processed the ADI 
algorithms with the same parameters. Note that for LOCI, we 
injected the fake planets in the cleaned and recentered datacubes 
before applying the reduction. Then the comparison between the 
injected flux to the retrieved one on the final fake planets im- 
ages was done by aperture photometry. This allowed to derive 
the actual attenuation for all separations from the central star by 
interpolating between the points. Finally, the 5cr detection limits 
were derived by taking the flux loss and the transmission of the 
neutral-density filter into account, and were normalized by the 
unsaturated PSF flux. 2D contrast maps were therefore available 
for each star, with each reduction techniques. 
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Table 6. Relative astrometry and photometry of the new visual 
binaries resolved with VLT/NaCo L' and ADI imaging mode. 



Name 


Date 


Sep. 


PA 


A L' 






(arcsec) 


(deg) 


(mag) 


HIP9685B 


11/20/2011 


0.303 ±0.013 


242.26 ± 1.16 


2.7±0.1 


HIP38160B 


1 1/25/2009 


0.141 ±0.013 


117.08 ±2.28 


3.1±0.3 


HIP53524B 


01/11/2012 


4.540 ± 0.009 


319.04 ±0.7 


6.2±0.1 


HIP59315B 


07/27/2010 


0.36 ±0.01 


259.2 ± 2.5 


5.1±0.3 


HIP88399B 


07/29/2011 


6.439 ± 0.009 


88.8 ± 0.7 


4.9±0.1 


HIP9358()B 


07/29/2012 


0.242 ± 0.013 


94.4 ± 1.3 


3.9±0.3 


HIP117452B 


07/12/2012 


3.667 ± 0.009 


237.8 ±0.8 


3.8±0.1 


HIP117452C 


07/12/2012 


3.402 ± 0.009 


238.6 ±0.8 


4.0±0.1 



To compare the detection performances between the stars, 
we built ID contrast curves. An azimuthaly averaging within 
1 -pixel annuli of increasing radius on the noise map was per- 
formed, followed by flux loss coiTection, and unsaturated PSF 
flux scaling. This approach however tends to degrade the perfor- 
mances at close-in separations due to asymmetric speckle and 
spider residuals on NaCo data, or even the presence of bright bi- 
nary component. To retrieve the detection performances within 
the entire FoV, we created composite maps between ADI pro- 
cessed and nADI processed ones. Indeed, beyond 2 " from the 
central star where the limitations are due to photon and read-out 
noises, nADI remains the most adapted reduction technique. It 
has been shown that the limiting long-lived (from few minutes 
to hours) q uasi-static speckles are well coiTelated for long-time 
exposures ("Maroi s et al.l 1200 6^ thus leading to a non Gaussian 
speckle noise in the residual image. Therefore, the definition and 
the estimation of cr to provide a detection threshold in the region 
limited by the quasi-static speckle noise might not be well ap- 
propriate and overestimated. However, the Gaussian noise dis- 
tribution being achieved in the background noise regime, the 5<j 
detection threshold corresponds to the expected confidence level. 
Moreover, the conversion from contrast to mass detection limits 
is much more affected by the uncertainties on the age of the tar- 
get stars and the use of evolutionary models than by uncertainties 
on the detection threshold. 



4. Results 

Our survey aims at detecting close-in young and warm giant 
planets, even interacting with circumstellar disks in the case of 
stars with IR excess. Four stars in the sample have been identi- 
fied as hosting substellar companions in previous surveys. The 
redetection of these companions allowed us to validate our ob- 
serving strategy, data reduction, and might give additional data 
points for orbital monitoring. We also imaged a transitional disk 
at an unprecedented resolution at 3.8 jum for the first time around 
HD 142527 (HIP 78092) which was presented in a dedicated pa- 
per (IRameau et al ] l2012h . However, we did not detect any new 
substellar companions in this study. 

In this section, we describe the properties of the newly re- 
solved visual binaries, we review the observed and characterized 
properties of known substellar companions and of the candidates 
identified as background sources. We then repoit the detection 
performances of this survey in terms of planetary masses ex- 
plored. Finally, we briefly summarize the results on some previ- 
ously resolved disks, especially about HD 142527, in the context 
of a deep search for giant planets in its close environment. 



4.1. Binaries 

Despite the rejection of known binaries with 1. - 2. " separa- 
tion, 8 visual multiple systems were resolved (Figure |3). Their 
relative position and magnitude are reported on Table |6] 4 pairs 
are very close-in, with separations below 0.4 " whereas the re- 
maining ones lie in the range 4-7". Only HIP 38160 has been 
observed at a second epoch and confirms as a comoving pair 
HIP 88399 B and HIP 1 17452 B and C were known from liter- 
ature and HIP 59315 B might be indeed bound to its host-star 
based on archive data. 

HIP 9685 - HIP 9685 is referenced as a A yu astrometric bi- 
nary (Makarov & Kaplan 2005) and was associated to a ROSAT 
source by IHaakonsen & Rutledgel (120091) . In this work, we re- 
port the detection of a close-in binary candidate companion at a 
projected separation of 14 AU if we adopt 47.7 pc of distance. In 
the 2MASS images in JHK taken in October 1999, a point source 
is visible toward the Noith-East direction, at a separation around 
12 " and a position angle of ^ 15 deg. From the two relative 
positions, it came out the 12"-CC in the 2MASS images is not 
compatible with a background star at the position of our 0.3 "- 
CC. It is also very unlikely that the 2MASS 12"-CC has travels 
in projection from 540 AU to 14 AU in ten years. Instead, the 
existence of the astrometric acceleration suggests that our NaCo 
0.3 "-CC is bound and is responsible for the astrometric signa- 
ture. The 2MASS PSF being symmetric, our 0.3 "-CC may have 
been at much smaller separation in 1999 since the orbital motion 
is significant, which does not contradict the proposed status. If 
it is true, we derive the mass of our 0.3 "- CC from the meas ured 
Ml' =4.1 mag using the isochrones from'Siess et al. V2OOO"), as- 
suming a solar metallicity, and an age of 30 Myr, to be 0.8 Mq, 
matched with a K6 star. 

HIP 38160 - A companion with a magnitude of L' - 7.84 
is present at two different epochs (2009-november and 2010- 
december) 4.8 AU away (0.141 ") from HIP 38160 at 35 pc. 
The companion shows a common proper motion with the central 
star between the two epochs. The 2MASS JHK images taken 
in 2000 also reveal an asymetric PSF which tends to confirm 
the bound status of this companion. According to ISiess et all 
(2000) isochrones for pre- and main-sequence stars, this com- 
panion should be of 0.6 - 0.7 Mq, assuming an age of 200 Myr 
and a solar-metallicity. Hence, HIP 38 160 B could be a late K or 
ealy M sta r HIP 38160 was already catalogued as an astromet- 
ric binary (iMakarov & K aplan 2005) and as a double-star sys- 
tem in t he Catalog of Component of Double or Multiple stars 
(CCDM iDommanget & Nv3[2003) . However, with a separation 
of 23.3 arcsec, this additional candida te turns to be only a visual 
companion (WDS. lMason et al.ll200lh . 

HIP 53524 - HIP 53524 lies at a very low galactic latitude 
{b < 10 deg). It is therefore very likely that the candidate com- 
panion, located at a large separation from the central-star (> 4"), 
is a background star Indeed, from //ST/NICMOS archive data 
taken in 2007, we measured the relative position of the well seen 
CC. Even not considering the systematics between the two in- 
struments, the CC turns out to be a background object. 

HIP 59315 - The star HIP 59315 is not catalogued as being 
part of a multiple physical system. However, it lies at relative 
low galactic latitude (Z? = 13 deg) but only one point source 
has b een detected w i th VU T/NaCo in ADI and L' imaging in 
2010. IChauvin et all (1201 Oh observed this star with VLT/NaCo 
in H band in coronographic mode and have identified an addi- 
tional background source more than 5.5 "away with a PA around 
100 deg. If our 0.36"-CC is a background contaminant, it would 
lie in April 2004 at a separation of 1 .35 "and a position angle of 
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Table 7. Relative astrometry and photometry of the known sub- 
stellar companions observed with VLT/NaCo L' and ADl imag- 
ing mode. 



HIP 9635 



HIP 38180 




Name 


Date 


Sep. 


PA 


AL ' 






(arc sec) 


(deg) 


(mag) 


HR 7329 b 


08/13/2011 


4.170 ±0.009 


167.43 ±0.7 


6.7 ±0.1 


AB Pictoris b 


1 1/26/2009 


5.420 ± 0.009 


175.2 ±0.7 


7.0 ±0.1 


jS Pictoris b 


09/27/2010 


0.383 ±0.11 


210.28 ± 1.73 


7.8 ±0.3 


HR 8799 b 


08/07/2011 


1.720 ±0.025 


62.9 ± 1.3 


9.9 ±0.1 


HR 8799 c 


08/07/2011 


0.940 ± 0.016 


321.1 ± 1.5 


7.1 ±0.2 


HR 8799 d 


08/07/2011 


0.649 ± 0.016 


207.5 ± 1.4 


7.1 ±0.2 



from the 2MASS survey. Given the separation in our observation 
and the L' magnitude, the CC is indeed the M dwarf companion, 
lying at 310 AU from the primary. 

HIP 93580 - The star is 70 Myr old A4V star at 55.19 
pc and b - -68.6 deg. A point source 3.9 magnitudes 
fainter than the primary is detected at a projected separation of 
13.24 AU. Neither archive nor second-epoch observation s could 
infer the status of this CC. Comparison to the lSiess et al. 
isochrones at 70 Myr with a solar metallicity would place this 
object as beeing an early M-type dwarf, with a mass of 0.5 Mq. 

HIP 117452 - Already known as a triple system 
(iDe Rosa et alj 1201 Ih from observations taken in 2009, the 
two companions are detected from our data in july, 2012. The 
brightest companion lies at ^ 160 AU in projection from the 
primary while the third component is at a separation of 1 1 AU 
from HIP 117452 B . The large error bars on the astrometry in 
iDe Rosa et al.l (|201 li) make difficult to infer any orbital motion 
of both companion in two years. 



Fig. 3. Resolved visual binaries with VLT/NaCo in ADI-L' 
imaging mode. HIP 9685, HIP 53524 and HIP 93580 have not 
been confirmed as double systems from second-epoch observa- 
tions or archived data. North is up, east to the left. Counts are 
displayed in linear scale but differently for each panel. 



254 deg so that it would have been detected on NaCo H images. 
The other possibility is that our 0.36"-CC is indeed bound to HIP 
59315 and was occulted by the mask. Therefore, it is likely that 
it is indeed bound to the star. This would imply for the compan- 
ion a projected separation of 14 AU and an absolute magnitude 
Mjf = 8.2 at 37.8 p c. The mass derived from the COND model 
dBarafFe et al.ll2003 ) assuming an age of 100 Myr is O.IM©, con- 
sistent with a late M dwarf. 

HIP 88399 - HIP 88399 is referenced as a double star in 
SIMBAD with a M2 star companion (HIP 88399 B) at 6.35 " 



We also have two spectroscop ic binaries (HI P 101800 and 
HIP 25486) in our survey. Pourbaix et al.l (12004 ) give a period 
of about 1 1 d, an eccentricity of 0.23, and a velocity amplitude 
of the primary of_Kj = 26 km/s for HIP 101800. For HIP25486, 
[kolmberg et ak] (| 200? |) report a standard deviation for the RV 
signal of 4.4 km/s, a SB2 nature with an estimated mass ratio 
of 0.715. However, we do not detect any source with a contrast 
from 5 mag at 100 mas up to 12 mag farther out of 1 .5" around 
both stars. Both companions are likely too close to their primary 
for being resolved, or even behind them. 

4.2. Substellar companions 

Four targets in the sample - HR 7329, AB Pictoris, /3 Pictoris, and 
HR 8799 - have previously reported to host a brown dwarf and/or 
planet companions ( Lowrance et al. 2000; Chauvin et al. 20051 
iLagrange et"aLll2010t iMarois et al. 2008 U2010I) . Only one iden- 
tified substellar CC to HIP 79881 has also been stated as back- 
ground object. We review below the latest results about these 
companions since their initial confirmation. Table |7] lists their 
relative astrometry and photometry from our observations ( see 
Figure nil. 

HIP 79881 - Clearly seen in 2010 july observations in L', 
the 11.6 mag-contrast CC to HIP 79881 (separation of 4.528 + 
0.008" and a position angle of 175.54 + 0.8 deg) has been also 
resolved in Keck/NIRC2 images in 2003 and 2005. The relative 
positions of the CC monitored for 7 years clearly showed that it 
is a background object. 
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HR 7329 AB Pic 




Fig. 4. Known substellar companions of HR 7329, AB Pic (Top 
left and right), and /? Pictoris and HR 8799 (Bottom left and 
right) observed with VLT/NaCo in ADI-L' modes. North is up, 
east to the left the count scale is linear. Note that, we could not 
retrieve with a good SNR HR 8799 e due to the low parallactic 
angle amplitude (20.9 deg). 



HIP 95261 I HR 7 3 29 - HR 7329 b was discov- 
ered by iLowrance et al.l (|2000|) with the Hubble Space 
Je/escope/NICMOS. It is separated from its host star, a mem- 
ber of the P Pic, moving group which harbors a debris disk 
dSmithetalJlIOOl) . of 4.17 ± 0.09 " (^ 200 AU at 48 pc) a 
position angle of 167.4 + 0.7 deg. The age and the known dis- 
tance of the star together with HSTISTIS spectra and photome- 
try from H to L' bands are consistent to infer HR 7329 B as a 
young M7-8 brown dw aif with a mass between 20 and 50 My. 
iNeuhauser et al.l ( 1201 lb conduct a 11 yr followed-up to confirm 
the status of the companion and try to constraint the orbital prop- 
erties. Due to the very small orbital motion, they concluded that 
HR 7329 B relies near the apastron of a very inclined - but not 
edge-one - and eccentric orbit. Our observations are consistent 
with the previous ones and exclude the presence of additional 
companions down to 4 My beyond 40 AU. 

HIP 30034/ AB Pic - This member of the Columba asso- 
ciation hosts a companion at the planet-brown dwa rf boundary 
of 13 + 2 My, discovered bv IChauvin et all (l2005h . Located at 
5.4 + 0.07" and 175.2 + 0.7 deg, AB Pic b has a mass of 13 - 
14 My deduced from evolution ary models and JHK photome- 
tries. Later on, iBonnefov et al.l (2010) and recently Bonnefoy 
et al. (2012, submit.) conduct observations with the integral 
field spectrograph VLT/SINFONl to extract medium-resolution 
{R,\ - 1500 - 2000) spectra over the range 1.1 - 2.5 nm. They 
derive a spectral type of LO-Ll, an effective temperature of 
- 1700 - 1800K, a surface gravity of log(g) ^ 4.5 dex by com- 
parison with synthetic spectra. The relative astrometry and pho- 
tometry we measured from our observations are similar to the 
previously reported ones. Further investigations are mandatory 
to derive, if similar, similar conclusions as for HR 7329. Due 
to our highest spatial resolution and sensitivity, surely planets 
more massive than 3 My can be excluded with a semi-major axis 
greater than 80 AU. 



HIP 27321 / p Pictoris -/3 Pic b (iLagrange et al.ll2010h re- 
mains up to now the most promising case of imaged planet prob- 
ably formed by core accretion. Recent results by (Chauvin et aT| 
l2012h . including measurements from this survey, refined the or- 
bital parameters with a semi-major axis of 8 -9 AU and an eccen- 
tricity lower than e < 0.17. In addition. iLagrange et al.l (1201 2ah 
could accurately show that the planet is located into the second- 
warped component of the debris d isk surrounding the star, which 
confi rms previous studies (Moui llet et al.|[T997t lAugereau et al.l 
12001 ) suggesting that the planet pl ays a key role i n the m or- 
phology of the disk. More recently, iLagrange et al.l (1201 2bb di- 
rectly constrain the mass of the planet through eight years high- 
precision radial velocity data, offering thus rare perspective for 
the calibration of mas s-luminosity re l ation of young massive 
giant planets. Finally, Bonnefoy et al.l ('2013') build for the first 
time the infrared spectral energy distribution of the planet. They 
derive temperature (1600 - 1800 K), log g (3.5 - 4.5), and lumi- 
nosity {log(L/LQ - -3.87 + 0.08) for/3 Pic b from the set of new 
and already published photometric measurements. They also de- 
rive its mass (6 - 15.5Mj) combining predictions from the latest 
evolutionary models ("warm-start", "hot-start") and dynamical 
constraints. 

HIP 114189 I HR 8799 - HR 8799 is a well-known A 
Boo, y Dor star, surrounded by a debris disk dPatience et al.l 
'201 1) and belonging to the 30 Myr-old Columba association 
(Zuckerman et al. 2011). It ho sts four planetary-mass compan- 
ions between 14 and 68 AU (iMarois etal.ll2008l 1201 Ol) which 
awards this multiple planet system being the first imaged so 
far. Spectra and p hotometry studies (e.g. iBowler et al.l 1201 Ol : 
'Janson et al. 2010) inferred those planets to rely between 5 and 
7 My. Soummeret al. (2011) monitor the motion of the planets 
b, c and d thanks to HST/NICMOS archive giving 10 yr am- 
plitude to constrain the orbits of these planets. Invoking mean- 
motion resonances and other assumptions for the outer plan- 
ets, they derive the inclination of the system to ~ 28 deg. 
Esposito et al.l (1201 2|) consider also the planet e for the dy- 
namical analysis of the system. They show that the copla- 
nar and circular system cannot be dynamically stable with the 
adopted planet masses, but can be consistent when they are about 
2 Mj lighter. In our images, HR 8799 b, c and d are clearly de- 
tected. The measured contrasts between the host star and each 
planet are very similar with the previously reported ones. No 
new orbital motion for planet b, d, and d is found from our ob- 
servations compared to the latest reported astrometric measure- 
ments by then end of 2010. Finally, the e component could not 
be retrieved with high signal to noise ratio due to the short am- 
plitude of parallactic angle excursion (20.9 deg). 

4.3. Detection performances 

Typical contrasts reached by our survey using ADI algorithms 
and nADI algorithms are presented on Figure |5] Note that us- 
ing cADI/sADI/rADI or LOCI, the performances are very sim- 
ilar, except within the exclusion area of LOCI. The median az- 
imuthally averaged L contrast vs the angular separation is plot- 
ted together with the best curves. The typical range of detection 
performances at all separations beyond 0.5 "is about 2 mag with 
a median contrast of 10 mag atO.5", 11.5 mag at l",and slightly 
below 12 mag at 2". Our best performances even reach very 
deep contrast, up to 13. 5 "at 0.5 "around HIP 118121. 

The detection limits (2D-maps and ID-curves) were con- 
verted to absolute L magnitudes using the target properties and 
to predicted masses using the COND03 (iBaraffe et all l2003l) 
evolutionary models for the NaCo passbands. 
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Semi-major axis (AU) Semi-major axis (AU) 



Fig. 6. Mean detection probability map of our survey as a function of the mass and semi-major axis (left) and probability curves for 
different mass bin as function of the semi-major axis (right). The mean is obtained over all targets of the survey. The detectability 
of the simulated planets was compared to the detection maps from LOCI and nADl algorithms and using COND03 (IBaraffe et al.l 
l2003h evolutionary models. 
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Fig. 5. Survey detection limits in L'-band contrast between the 
central star and any point source vs the angular separation, at 
the 5cr level, using VLT/NaCo in ADl mode. Solid lines are rep- 
resentative of median performances whereas the bottom dashed 
ones are for the best performance reached by our survey. cADl 
(red) and LOCI (blue) are very similar whereas sADI (orange) 
and rADI (aqua) remain slightly above. Note that the LOCI 
curves stop at 0.15 "due to the central exclusion area. 



The overall sensitivity of our survey can be estimated us- 
ing Monte-Carlo simulations. We use an optimized version of 
the MESS code (Bonavita et al. 2012) to generate large popula- 
tions of planets with random physical and orbital parameters and 
check their detectability by comparing with the deep detection 
limits of our survey. We performed simulations with a uniform 
grid of mass and semi-major axis in the interval [1,20] Mj and 
[1, 1000] AU with a sampHng of 0.5 Mj and 2 AU. For each 
point in the grid, lO** orbits were generated. These orbits are ran- 
domly oriented in space from uniform distributions in sin(/), a>, 



Q., e < 0.8, and rj^. The on-sky projected position (separation 
and position angle) at the time of the observation is then com- 
puted for each orbit. Using 2D informations, one can take into 
account projection effects and constrain the semi-major instead 
of the projected separation of the companion. We ran these sim- 
ulations for each target to compute a completeness map with no 
a piioii information on the companion population and therefore 
considering a uniform distribution in mass and semi-major axis. 
In this case, the mean detection probability map of the survey 
is derived by averaging the 59 individual maps. The result is il- 
lustrated on Figure |6] with contour lines as function of sma and 
masses. Note that the decreasing detection probability for very 
large semi-major axis reflects the fact that such companion could 
be observed within the FoV only a fraction of their orbits given 
favorable parameters. The peak of sensitivity of our survey oc- 
curs at semi-major axis between 40 and 300 AU, with the high- 
est sensitivity around 100 AU. Our survey's completness peaks 
at 94 %, 80 %, and 58 % for 10 Mj and 5 Mj at 100 AU, and 
3 Mj at 220 AU, respectively. The overall survey's sensitivity at 
1 Mjis very low, with a maximum of 15 % at 141 AU. 

4.4. Giant planets around resolved disks 

Among the targets of our sample, HD 142527 remains the 
youn gest one. It has been identified bv lAcke & van den Anckeil 
(12004) as member of the very young (^ 5 Myr) Upper 
Centaurus-Lupus association. HD 142527 was selected to take 
advantage of the capability offered by VLT/NaCo in thermal 
and angular differential imaging to resolve for the first time at 
the sub arcsecond leve l its circumstellar environment. Indeed, 
iFukapawa et alj (l2006h detected a complex transitional disk in 
NI R with a huge g ap up to 100 AU. Our observations reported 
by i Rameau et al.l (12012 ') confirm some of the previously de- 
scribed structures but reveal important asymmetries such as sev- 
eral spiral arms and a non-circular large disk cavity down to at 
least 30 AU. The achieved detection performances enable us to 
exclude the presence of brown dwarfs and massive giant plan- 



^ They correspond to respectively the inclination, the argument of 
the periastron with respect to the line of nodes, the longitude of the 
ascending node, the eccentricity, and the time of passage at periastron. 
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Fig. 7. Detection limits of giant planets in Jupiter mass vs the 
projected distance from the central star for stars with known 
debris disks. The detections are expressed at the 5cr level, us- 
ing VLT/NaCo in ADI mode, after LOCI (close-in region) and 
nADI (background limited region) processing. The COND03 
dBarafFe et al.l2003 ) mass-luminosity relationship is used to con- 
vert from L' contrast to Jupiter mass. Note that for AU Mic, the 
used model does not go under 0.5 Mj, so we cut the curve when 
it reaches this limit. 

ets beyond 50 AU. In addition, two sources were detected in 
the FoV of NaCo. The relative astrometry of each CC was com- 
pared to the one extracted from archived observations and to the 
track of the relative position of stationary background object. 
Both were identified as background sources. Finally, the candi- 
da te substellar compa nion of 0.2 at 13 AU recently reported 
bv lBiller et all (12012') may be responsible for the structured fea- 
tures within the disk. These structures might also be created by 
type II migration or by planetary formation through GI. Further 
investigations of the environment of HD 142527 could set con- 
straints on planetary formation and disk evolution. 

In our sample, we also have included some stars surrounded 
by debris disks that have been resolved in previous observations 
(49 Cet, HD 10647, HD 15115, Zeta Lep, 30Mon,HD 181296, 
HD 181327, HD 181869, HD 191089, and AU Mic). Both obser- 
vation and reduction processes having been designed to search 
for point sources, we do not report results about disk properties. 
However, we investigate the presence of giant planets and plot 
in the Figure |7] the detection limits about these stars in terms 
of mass vs projected separation. These limits could help to con- 
strain some disk properties which can be created by gravitational 
perturbation of giant planets. The detection sensitivity around 
AU Mic reaches the sub Jovian mass regime at very few AUs 
from the star because AU Mic is a very nearby M dwarf. On the 
other hand, the detection limits around HD 181869 are not very 
good due to bad quality data. We remind that these limits are az- 
imuthally average so that there might be affected by the presence 
of the disk. 

5. Giant planet properties, occurrence and 
formation mechanisms 

The frequency of giant planets / can be derived using known 
planets and detection limits in case of a null detection. For 
an arbitrary giant planet population, one can compute within 
the mass and semi-major axis ranges probed by the survey. 




10 20 30 40 50 

Frequency (%) 



Fig. 8. A-F sample probability density of the fraction / of stars 
hosting at least one giant planet at wide orbit, taking into ac- 
count the detections {{dj]) of the two planetary systems /3 Pictoris 
and HR 8799 and a linear-flat prior The interval considered is 
[1, 1000] AU and [1, 13] Mj. The confidence interval at 68 % 
confident level, labelled as CL, (blue) and at 95 % CL (dark blue) 
are over plotted to the distribution. A uniform planet population 
has been generated with random orbital parameters. 

Numerous deep imaging surveys did not report the detection of 
at le ast one substellar or planetary mass companion. The authors 
(e.g. Kasper et al. ''2007^, 'Laf reniere et"aLl l2007f: 'Niels en & Cioig 
2010; Chauvin et al. 2010) nevertheless performed statistical 
analysis with MC simulations to fully exploit the potential of 
thei r data and p r ovided upper limits to the frequency of plan- 
ets. IVigan et aP (1201 2h tooks into account the planets already 
identified (J3 Pic, HR 8799) to derive also lower limits to this 
frequency. 

In this section, we derive the rate of wide-orbit giant plan- 
ets following the statistical approach used in previous works 
( ICarsonetalJ r2006; Lafreniere et al. 2007; Vigan et al. 201^ 
and d escribed in the appendix. Similarly to Bonavi ta et all 
(1201 2h . we take into account the binary status of some stars to 
exclude semi-major axis values for orbits which would be un- 
stable. The whole section relies on the two previously defined 
sub-samples : 37 A-F stars and 29 A-F dusty stars (Table |4). 
First, the frequency of wide-orbits planets is derived assuming a 
uniform distribution. We then use and discuss the extrapolation 
of RV statistics to wide-orbits in light of DI planets. Planet for- 
mation is finally considered, GI and CA, to estimate the impact 
on the observed occuiTence of giant planets. 

5.1. Occurrence of giant planets from a uniform distribution 

Assuming a uniform distribution of planets in a giid 
[1, 1000] AU and [1, 13] Mj, we use MC simulations to measure 
the detection probability pj around each star given the detection 
sensitivity, as in section 1431 The probability density function is 
then derived using equations lA. 1 1 and l A.2I assuming a flat prior 
Finally, the confidence interval of the true / is computed using 
equation lA.4l 

First, we focus on the relevant A-F statistical sample (see 
Table[TJ. In this sample, 2 stars harbor at least one giant planet, 
/3 Pictoris and HR 8799. Since these two stars match our selec- 
tion criteria, they were originally included in our sample, even 
if the giant planets have been discovered by other observations. 
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We thus take, in our analysis, at least two planetary detections 
(since our observations lead to the confirmation of their status). 
The Figure|8]shows the posterior distribution as function of / in 
the interval [1, 1000] AU and [1, 13] Mj. The observed rate of 
giant giant planets at wide orbit leads / to be 16.1^j^ 2 % with a 
confidence level (CL) of 95 %. At 68 % CL, this rate becomes 
16.1^5 3 %. If one considers the sample of 29 A-F dusty stars, 
thus with the same planet detections, the giant planet occurrence 
is 21.4+^^-^ % at 95 % CL or 21 .4":)^ % at 68 % CL. Due to 
our poor sensitivity to close-in and/or low mass planets, these 
values are relatively high. If we re strain the int erval of interest 
to [5, 320] AU and [3, 14] Mj as inVigan etal] (12012 ). then the 
A-F sample has an occurrence of giant planets of 7.4;^2 4 at 
68 % CL which matches the results obtained by the authors. 

The same approach can be done to derive the frequency of 
brown dwarfs. Taking into account the detection of HR 7329 
b in the A-F sample, / is 6.5!^<;f % and 8.7!fg'* % in the A- 
F dusty sample in the interval [i, 1000] AU and [14,75] Mj at 
95 % CL. The confidence interval is smaller compared to the 
statistical results for giant planets due to our high sensitivity to 
brown dwarfs. 

Finally, the full survey of 59 young, nearby, and B- to M-type 
stars can also give some constraints on the occurrence of planets 
within a broad sample of stars. Since the companions to AB Pic, 
HR 7329, HR 8799, and /3 Pictoris were previously detected out 
of our observations, we consider here a null detection within 55 
stars of our survey. Using equation lA.6l an upper limit to the fre- 
quency of giant planets can be derived with our detection limits 
and MC simulations. It comes out that less than 25 % among 
our 55 stars harbors a giant planet in the range [40,600] AU 
and [5, 13] Mj at 95 % CL. Note that this upper limit shaiply 
increases towards smaller mass planets and also to a wider semi- 
major range due to our poor sensitivity. We also remind that this 
sample is statistically less relevant than the previous ones since 
it is more heterogenous in terms of stellar mass, distance, and 
spectral type. 

5.2. Giant planet population extrapolating radial velocity 
results 

Radial velocity results provided a lot of statistical results on the 
giant planet properties but also on the distribution of the popula- 
tion with respect to the mass and/or semi-major axis. However, 
such results are intrinsically limited so far to close-in planets 
(typically 3-5 AU). Numerous publications present statisti- 
cal analysis on giant planet detected by deep imaging using 
extrapolation of the RV frequencies and distri butions to plan- 
ets on larger orbits (j Eafreniere et al. 2007; Ka sper et al.ll2007t 
iNielsen & Closell2010riVigan et alJ i2012). We briefly present in 
the following sections the outcomes of our sample based on the 
same approach, considering the detections around /3 Pic and HR 
8799. 

5.2.1 . Extrapolation of the radial velocity planets distribution 

We assume that the mass and semi-major axis distributions 
follow the simple parametric laws of index a and /3 : dN oc 
M"„dM„ with - 1.31 and dN k a^^da with yS = -0.610 

(ICumming et alJl2008h . Here, we blindly extrapolate the distri- 



While they derived the distribution for mass and period in logarith- 
mic bins using a as the index for P, we used the mass and semi-major 
axis distribution with linear bins using a referring to a. 
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Fig. 9. Contours showing the confidence level at which we can 
reject a planet distribution following a mass power law dN oc 
M"dM and a semi-major axis one dN oc d^da at a given semi- 
major axis cutoff. Top : Semi-major axis cutoff vs semi-major 
axis power law index fi for a - -1.5. Bottom : Semi-major axis 
cutoff vs semi-major axis power law index fi for a - 1.1, i.e. with 
more massive planets. The upper cutoff for the ICumming et all 
( l2008h P index of-0.61 is over plotted (dotted line). These fig- 
ures use LOCI and nADI detection limits and the COND03 
jBaraffe et alj|2003h for the mass-luminosity relationship. 



bution to larger semi major axis while it is formally valid only 
for planets with semi-major axis below ^ 3 AU. 

For this calculation, we populate a grid of mass and semi- 
m ajor axis in the iri t ervals [1, 13] and [1, 1000] (Mj and AU) to 
the lCumming et all (l2008h power-laws and run MC simulations 
to derive the probability density distribution as in section |5] If 
the giant planet population on wide orbits follows the RV power- 
laws, then their frequency range from our study is 22.0!]^-^% at 
95 % CL or 22.0^_]\^ % at 68 % CL in the A-F sample. This rate 

becomes 28.3!^^-^ % at 95 % CL or 28.3+J'*/ % at 68 % CL in 
the A-F dusty sample. 

However, there are intrinsic limitations on this study and the 
outputs, eventhough close in values as the ones reported among 
an uniform distribution, have to be taken with care. The used dis- 
tribution fits the statistic for solar-type stars up to few AU (com- 
ing from RV surveys) and is arbitrarily extrapolated to large sep- 
arations. There is also no evidence that the few planets detected 
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so far at large orbit separations have similar properties and dis- 
tributions. 

5.2.2. Constraining tine parametric laws for tine giant planet 
distribution 

This likelihood approach answers the question : 'How consistent 
is a given giant planet population with our observing results?' 
Answering this question requires 1/ to know all giant planet pop- 
ulation parameters and 2/ to know the fraction of stars with giant 
planets according to this given distribution. For each star, we can 
then derive the number of expected detections given the detec- 
tion sensitivity and compare to our observations. Such compari- 
son allows us to constrain a given distribution of wide orbit giant 
planets. Likewise, a giant planet population in which 95 % of the 
predicted planets would have lead to detections can be consid- 
ered as strongly inconsistent with our survey. Finally, this study 
relies on the strong assumption that we know the frequency of 
giant planets in the range where our survey is sensitive to. 

In the following, we use a pop ulation of plane t s give n by 
power laws similar to the ones from iCumming et al.l (l2008b and 
we also add an additional parameter which is acutoff , the semi- 
major axis beyond which there are no planets. Our intervals of 
interest for the simulation are [1, 1000] AU and [1, 13] Mj, nor- 
maUzed with / = 1 0.5 % over the range fO .3, 10] Mj, [2, 2000] 
days in period from 'Gumming et al.l (l2008l) . /norm is thus set with 
the ratio of the integrated power laws for a pair {a, /3) over 
[1, 13] Mj and [1, 1000] AU and the same over the RV ranges. 

We explored a grid of a, /3, and Ocutoif with a sampling of 
0.2 for the power law indices and 20 AU for the cutoff to de- 
rive the expected number of planets for each combination of 
parameters over the A-F sample (similar results are obtained 
with the A-F dusty sample). We illustrate the confidence level 
at which we can reject each model in Figure |9] as a function of 
j3 and Ocutoff for the A-F sample for two values of a : -1.5 and 
1.1, values corresponding to the extrema of our grid and thus 
showing the trend of the rejections. All results (ours and pre- 
vious publications) are consistent with a decreasing number of 
giant planets (P < -0.61) w hile their mass increase. Considering 
the Cummin g et al.l (|2008|) distributions, a semi-major axis cut- 
off around 45' - 65 AU at 95 % CL is found. 

We remind that mixing power-laws derived from RV and gi- 
ant planets with possibly different formation processes and evo- 
lutions has to be considered with caution. 



where Q is the Toomre parameter, c, the sound speed, k the 
epicyclic frequency and S the gas surface density. The Toomre 
parameter is fulfilled at larger radius only when the local mass is 
high enough. Therefore, fulfilling the Toomre criteria leads to a 
given value Z which can be converted to mass and thus states a 
lower limit in the mass-sma diagram : 

r 

where H - Cj/Q is the disk scale height. The other parame- 
ter which drives the instability is the cooling time, which, 
if higher than a few local keplerian timescale Q i.e. at 
small separation, stabilizes the disk through turbulent dissipa- 
tion (Gammie 2001; Rafikov & Goldreich 2005) . It thus puts an 
upper boundary in the mass-sma diagram : 

Mf = I.(2nHf 

where 2nH is the wavelength of the most unstable mode. Such 
boundaries, being global and excluding long term evolution, as- 
sume planets formed in-situ with masses of the disk fragments. 




Semi-major axis (AU) 



5.3. Giant planet formation by gravitational instability 

Gravitational instability is a competitive scenario to form giant 
planets, specially at large separations. Such a process becomes 
more efficient within massive disks, i.e. around massive stars. 
Since our statistical sample contains A-F and/or dusty stars, i.e. 
massive stars, we were strongly tempted to test the predictions 
of GI models wit h our observing resu lts. We hence adopted the 
same approach as iJanson et al.l (1201 ll) . The reader is refered to 
iGammid (l200lh for a detail description. The ID current model 
of disk instability provides formation criteria, which if fulfilled, 
create an allowed formation space in the mass-sm a diagram. The 
first one is the well known Toomre parameter ('Toomrel ll98ll) 
which has to be low enough to allow local gravitational instabil- 
ity in a keplerian accretion disk : 



Fig. 11. Disk instability model predictions for /3 Pictoris. The 
lower solid curve corresponds to the Toomre criteria which ex- 
cludes the formation of planets below it. The upper dashed curve 
corresponds the cooling criteria which only allows disk fragmen- 
tation below it. The allowed formation space is in between. The 
ID average detection limit curve has been over plotted (blue 
solid line) but with a projected separation. The location of /3 
Pictoris b has been over plotted to the graph with new error bars 
from .Bonnefov et al.. (.2013.) . 

The model computes both boundary curves for each star in 
the sample, taking into account the stellar mass, luminosity, and 
metallicity, the later being extracted from lAmmons et aT l (l2006l) 
or set to the solar one when the information was not available 
and luminosities derived from isochrones of ISiess et al.l (l2000l) 
using their absolute K magnitude, spectral type age, and metal- 
licity. The model is very sensitive to the stellar luminosity since 
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Fig. 10. Left : Mean detection probability map of the A-F sample as function of the mass and semi-major axis of the substellar com- 
panion. A uniform grid distribution has been used to generate the population but the formation limits derived from the gravitational 
instability models exclude for each star the planets which do not fulfill both criteria. The detection maps from LOCI and nADI 
algorithms were used with COND03 (Baraffe et al. 2003) evolutionary models to convert from contrast to mass. Contour lines are 
regular from 0.5 to 0.9 plus one at 0.1. Right : Corresponding estimation of the upper limit, with a confidence level of 95 %, on 
the fraction of stars from the A-F sample, harboring at least one object companion in the same semi-major axis range. The curves 
for 25 Mj, 18 Mj, 13 Mj and 10 Mj have been plotted since lower mass planets are not allowed to form via the disk instability 
mechanism. 



strong illumination favors t he disk to be gravitationally stable 
jKenvon & Hartmannlll987l) . Figure [TT] shows one example for 
P Pictoris. The Toomre and cooling criteria are fulfilled around 
40 AU for very massive planets (> 3 Mj) and this trend rapidly 
increases with the separation, thus leading to the brown dwarf 
and stellar regimes. Note that considering a lower mass star 
would lead to push the boundaries inwards. 

We then run MC simulations in a uniform grid of mass and 
semi-major axis in the interval [1, 100] and [1, 1000] (Mj AU) as 
in section |5] Points of the grid out of the allowed range for each 
star are removed according to the formation limits. We remind 
that these predictions are not normalized due to the absence of 
knowledge on physical and statistical properties of protoplan- 
etary disks in which Gl starts. The mean detection probability 
over the A-F sample is plotted on Figure [TO] left panel. Only 
high mass planets and brown dwarfs fulfill the formation criteria 
and there are almost all detectable. 

/? Pictoris b, HR 8799 b, c, and d are too light and too close 
to their stars, so they do not fulfill both Gl boundary conditions 
(M < 10 Mj below 70 AU). Therefore, we cannot use these de- 
tections to derive the rate of giant planets according to Gl mech- 
anism. We instead estimate the upper limit on /, /max using equa- 
tion lA.6l We derive and plot /^ax (Figure[TO] right panel) for the 
A-F sample only for the mass regime allowed by this approach, 
which extends between very few tens of AU. The curves are off- 
set one from another due to the fact that higher mass object can 
be formed in-situ at larger distance from the central star. It came 
out that less than 20 % (< 23 % for the A-F dusty sample) stars 
harbor at least a 13 Mj planet between 40 and 60 AU and less 
than 25 % (< 27 %) a 10 Mj in the range [32,45] AU. 

On the other hand, 1/ Figure [TOl left panel, shows our high 
sensitivity to brown dwarf on wide orbits, and 2/ HR 7329, be- 
longing to the A-F dusty sample so as to the A-F one, hosts a 
detected brown dwarf companion for which the formation is al- 
lowed according to our Gl model. We can therefore estimate the 
rate of formed objects as in section |5] Since Gl can form plan- 
etary to brown dwarf mass objects, we explore the full range 



[1,75] and [1, 1000] (Mj and AU). We found that / equals to 
3.2+^-,2 % for the A-F sample and 4.3!^^ % for the A-F dusty 
one at 68 % CL if formed by this mechanism. 

It comes out that such Gl boundaries prevent the formation 
of low mass and close-in giant planets but would enhance the 
presence of brown dwarf and low mass star companions. Since 
high mass stars would facilitate the Gl mechanism by harboring 
massive disks, one would expect to find a higher occurxence of 
low mass stars or substellar companions rather than planets and 
a continuous distribution between the wide orbit giant pla nets 
detected so far and higher mass objects (iKratter et al.ll2010l) . 

This approach is a first step towards understand planet for- 
mation by Gl and the analysis can be improved by taki ng into 
account the following steps. First, iMeru &Batel (1201 Ih shows 
that using proper 3D global radiative transfer codes and hydro- 
dynamical simulations, closer-in disk region might become un- 
stable, phenomen a which was prevent ed assuming global simple 
cooling time law. lKratter & Murrav-Cl av (201 1 ) refines the def- 
inition of Q and the cooling time leading Gl to be possible at 
smaller separations. Second, the probability of clump formation 
towards planets was assumed to be one but long h ved clumps re- 
quire careful considerations about dis k dvnamics (| D urisen et aP 
l2007h . Then, clump evolution (e.g. lGalvagnietal.ll2012h and 
fragmentation might lead to the formation of lower mass plan- 
ets. HR 8799 seems a good test-case for such hypothesis. Indeed, 
the three outer planets orbit the star too far away to have form 
via core accretion. Gravitational instability naturally comes out 
as alternative scenario. However, each planet, with its mass and 
separation, does not fulfill the Toomre and cooling time crite- 
ria following our models. Considering all three together, even 
four mass planets 30 Mj) onto a single disk fragment at a 
mean separation satisfies our boundaries. One might speculate 
that this clump would have then broken after collapse leading 
to individual evolution of the planets. Finally, long term clump 
evolution was also not taken into account in our study. A self- 
graviting clump will still accrete a large amount of gas. Even 
if the disk fragment into an initially planetary mass clump, this 
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Fig. 12. Synthetic planetary population predicted assuming the 
core-accretion scenario similar as in Mordasini et al. (2012) for 
the case of 2 Mj central stars. The 20 %, 60 % and 80 % detec- 
tion probability curves (red) are overplotted to the planet popu- 
lation revealing the poor sensitivity of our survey on such kind 
of formed planets. The SPHERE/IFS performanc es (dashed) 
have also been plotted assuming contrast curves from lMesa et alj 
(1201 1 ). The short high-sensivity window is due to the small FoV 
of the IFS (1.77") which can be overcome with the larger IRDIS 
FoV. 



fragment will accrete gas, become more massiv e, and thus might 
exceed the deuterium burning limit mass (B0SS1I2OIII) . However 
such formation takes about 10^ yr, gas accretion is expected 
to be turned off by disk dissipation by strong UV irradiation 
of th e surrounding high m ass stars in the host-star forming re- 
gion dPurisen et al. I l2007t) so that one might expect light clump 
growth to stop before getting too massive. 

5.4. Giant planet formation by core-accretion 

W e now investigate the pl anet formation and evolutionary model 
of iMordasini et al.l (1201 2') which predicts the final state of plan- 
ets following the core accretion scenario and normalized by the 
frequency of observed disks. The synthetic population is calcu- 
lated assuming a 2 central star, a mean disk lifetime of 4 Myr, 
that gap formation does not reduce gas accretioifl and consid- 
ering one embryo per disk-simulation (hence no outward migra- 
tion of resonant pairs or scattering possible). A comparison with 
RV data shows that this simulation produces too massive and 
too close-in giant planets, but this synthetic popu lation can be 
considerer as a rough approximation (see also in lAlibert et al.l 
l201lh . We then try to test the predicted expected population at 
wide orbits with this approach so with direct imaging results. 
We run the MC simulations ran with planets extracted from this 
synthetic population. 10"* random orbits were generated for each 
planet and the projected position on the sky was computed as 
before. 

In Figure [12] we show the extracted planet population (al- 
ready normalized) as well as the detection probability at 20 %, 
60 %, and 80 % for the A-F sample. It comes out that there 
is no incompatibility between the synthetic population and the 



' This quest ion is still debated since gap formation might lead to a 
reduction (e.g. iLubow et al.lll999t) , but this might be coun terbalanced 
by the effects of eccentric instabiUty dKlev & Dirksenll2006h . 



results of our survey. Indeed, we are marginally sensitive to 
the farthest predicted giant planets (the predicted fraction with 
detectable planets is around 0.06 %) so we cannot reject their 
existence. Notwithstanding, CA is expected to become ineffi- 
cient to form planets at separations larger than a tens of AU. 
Moreover, the domain probed by our detection sensitivity, i.e. 
beyond 40 AU, well matched the region where CA inoper- 
ates as seen in Figure [12] The 5-20 AU gap between deep 
imaging surveys and those from radial velocity will be at least 
partly fill in thanks to the forthcoming extreme adaptive optic 
instruments VLT/ SPHERE (Beuzit et al. 2008) and Gemini/GPI 
([Macintosh et al. 2008) thanks to excellent detection limits and 
lower inner working angles. Using the same MC simulations, we 
compute the mean detection probability curves (Figure [I2]0 We 
show that the improved capabilities of SPHERE will allow in- 
deed to decrease this gap, by detecting a few Jupiter-like planets 
down to 5 AU. However, its overall sensibility (in mass and sepa- 
ration) will not allow to entirely probe the predicted giant planet 
population (the predicted fraction with detectable planets being 
around 0.6 %). Another complementary way to fill this gap is 
to use both RV and direct imaging on selected young targets, as 
demonstrated in Lagrange et al. (2012, subm.). 



6. Concluding remarks 

Here, we have reported the observations and analysis of a sur- 
vey of 59 stars with VLT/NaCo at L'-band (3.8/im) with the goal 
to detect and characterize giant planets on wide-orbits. The se- 
lected sample favors young, i.e. < 70 Myr, nearby, < 100 pc, 
dusty, and early-type stars to maximize the range of mass and 
separation over which the observations are sensitive. The opti- 
mized observation strategy with the angular differential imaging 
in thermal-band and a dedicated data reduction using various al- 
gorithms allow us to reach a contrast between the central star 
and an off-axis point source of 12 mag at 0.3", 13.5 mag at 0.5" 
up to 14 mag farther away in the best case. Despite the good 
sensivity of our survey, we do not detect any new giant planet. 
New visual binaries have been resolved, HIP 38160 confirmed 
as a comoving pair and HIP 79881 and HIP 53524 conffimed 
as background objects. We also report the observations of a per- 
fect laboraty-case! for disk evolution with the sub-arcsecond re- 
solved disk surrounding HD 142527 (dedicated publication in 
Rameau et al. (2012). 

We used Monte-Carlo simulations to estimate the sensitiv- 
ity survey performance in terms of planetary mass and semi- 
major axis. The best detection probability matches the range 
40 - 300 AU, with maxima at 93 % for a 10 Mj planet and 58 % 
for a 3 Mj planet. Brown dwarfs would have been detected with 
more than 70 % probability within the same semi-major axis 
range. 

A dedicated statistical analysis was carried out to under- 
stand and constrain the formation mechanism of giant planets. 
From literature and archive data, we focused on two volume- 
limited samples, representatives of almost 60 % to more than 
70 % of the full set of stars being younger than 100 Myr, 
closer than 65 pc, to the South (dec < 25 deg), A-or F-type, 
and with/without infrared excess at 24 and/or lOfim. We com- 
puted the frequency of giant planets at wide orbits, in the interval 
[1, 13] and [1, 1000] (Mj and AU), summarized in Table[8]: 



Only a small semi-major axis range is covered by the curves since 
we considered only the IFS instrument which has a small FoV (1.77"). 
Larger FoV will be provided by the IRDIS focal instrument. 
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Table 8. Confidence interval of the frequency of giant planets 
with a confidence level of 68 % reported in this work around 
young, nearby and dusty A-F stars assuming different planet 
populations, the detections of (3 Pictoris b and the system around 
HR 8799 for planets, and the detection of HR 7329 b as brown 
dwarf. 



Sep. range 
(AU) 


Mass range Frequency 
(Mj) (%) 


Distribution 


A-F sample 


[1,1000] 
[1, 1000] 
[1, 1000] 


[1,13] 10.8-24.8 
[1,13] 14.8-36.8 
[1,75] 2.2-5.4 


flat 

Cu08 

flat-^GI 


A-F dusty sample 


[1, 1000] 
[1, 1000] 
[1,1000] 


[1,13] 14.3-35.2 
[1,13] 18.7-47.9 
[1,75] 2.9-6.7 


flat 

Cu08 

flat+GI 



Notes. Results on the frequency of giant planets are reported accord- 
ing to a flat/uniform giant planet distribution or to a po wer law distribu- 
tion of giant planets with the mass and semi-major axis dCumming et alJ 
[2008) or driven by formation boundaries according to the gravitational 
instability scenario. The first two scenarii consider the detections of two 
planetary systems whereas the later set with one brown dwarf detection 
since planet formation through GI remains very low because they are 
not allowed to form closer-in and become too massive farther out. 

- in both A-F and A-F dusty samples, two giant planetary sys- 
tems have been detected so far : p Pictoris and HR 8799, 
yielding a wide-orbit giant planet occurrence between 4.9 % 
and 42.4 % for the A-F sample and between 6.5 % and 
56.9 % for the A-F dusty sample at 95 % CL, assuming a 
uniform distribution. These results are consistent with the 
upper limit found in the litterature and also with the rate of 
planets jiround the volume-limited sample of 42 A-type stars 
bv lVigan et alJ (120121) (see results in Table[T]i. 

- if the population of giant planets on wide orbits follows the 
distribution of the one s detected by RV below 5 AU from 
ICumming et al.l (l2008l) . the 95 % confidence interval for / 
is 6.7 - 59.5 % for the A-F sample, and 8.6 - 66.2 % for 
the A-F dusty sample. We remind that such an assumption 
is probably incorrect as it implicitly assumes that wide orbit 
giant planets have similar origins and properties as close-in 
ones around Sun-like stars. 

- planets formed via gravitational instability within protoplan- 
etary disks are expected to be massive and to orbit far away 
from their host stars if they remain in situ. We consider 
such planets to form and remain in situ where such insta- 
bility could occur Since /? Pic b, HR 8799 b,c, and d are 
not allowed to form via GI according to our model, we only 
computed upper limit to the frequency of giant planets. We 
find that less than 25 % of stars could form and retain a 
10 - 13 Mj between 30 - 60 AU in the A-F sample (so as 
for the A-F dusty sample). Closer in, disk instabilities are 
quickly prevented so that no planet can be formed, whereas 
the disk fragmentation rapidly leads to brown dwarfs and 
stellar regimes farther away. 

These results would coiToborate a correlation between the 
presence of debris disk and giant planets since the rates tend to 
be slightly eventhough higher in the A-F dusty sample than in 
the A-F one. They also point towards a similar occurrence of 



giant planets on wide (from AO imaging) and close (from RV 
measurement) separations. They suggest a break up of the pos- 
itive correlation between the separation, the mass, and the dis- 
tribution derived from close-in CA planets and the population of 
wide orbit giant planets. The later is consistent with a decreasing 
distribution with larger semi-major axis. All previous surveys re- 
ported such bimodal behavior of the distribution which would be 
a sign post of different modes of gas giant formation (Boley et al] 
|2009|) . Upcoming extreme AO surveys will probe the transition 
region between the two regimes revealing if it is co ntinuous, i.e. 
same formation process or, not (iKratter et al.ll20ldl) . 

Since our survey is very sensitive to high mass objects (i.e 
> 10 Mj), we can derive the rate of brown dwaifs to be 6.2^j ^ % 
in the A-F sample and 8.7^2^ % in the A-F dusty one, at 68 % 
CL in the interval [14,75] and [1, 1000] (IVlj AU) assuming a 
uniform distribution, results which are consistent with the lit- 
erature. From the GI formation boundaries, these rate become 
3.2+^-^ % in the A-F sample and 4.2+^-^ % in the A-F dusty one. 

We finally recall that all the detection limit estimations are 
based on mass-luminosity relations that are still debated. Also, 
they strongly rely on age estimates which are much less accu- 
rate for early-type stars when they do not belong to a moving 
group. The long term dyna mical evolutio n of plane tary system 
(migr ation in e.g. |P apaloizo u et al.ll2007l scattering ICrida et"al1 
2009; Ravmo nd et al.ll2012l) plays also a key role, and the giant 
planet distribution at a given age could be different from the one 
at formation stages. 
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Appendix A: Statistical formalism 

Ou r likelihood analys is approach follows the work done 
bv ICarson et alJ (l2006h : iLafreniere etaH (l2007l) : IVigan et alJ 



(I2012h . We nevertheless recall here the steps. 

The principle of detecting a planet around a star is a 
Bernoulli event. We note, pj the probability to detect a giant 
planet around a star j if it is indeed here. This probability de- 
pends on the distance, the luminosity, and the age of the star, and 
on the projected position, and the luminosity of a planet, and also 
on the instrumental performances. We also note the fraction of 
stars / harboring at least a planet in the interval [m^in, '«max] and 
[flmin, Omax]- We assume / to be constant around the star sample. 
For a given /' star, the probability of detecting a giant planet com- 
panion is fpj. Given the observational results of a survey, one 
can let dj to represent the detection efficiency such that dj - 1 
if a planet has been detecting around the star j and otherwise. 
Therefore, the likelihood function of the data for a set of star 
will be the product of each Bernoulli event since there are inde- 
pendent, so that : 

N 

mdj]\f) = n^-^^^^"'^^ - ^Pj'^''"' 1^ 

Then, we can apply the Bayes'rule which links the likeli- 
hood function of the data ({dj}) given the model / X,({dj]\ f) to 
the probability density of the model given the data, or posterior 
distribution P(f\{dj}). We get : 

mdi]\f)P(f) 

P(f\{dj}) = . ' ' (A.2) 
X £ildj]\f)P(f)df 

The Bayes' rule also remains on the assumption on the ini- 
tial probability of the model, or prior distribution P(f), which 
can be the most controversial part. One can use the posterior 
distribution from previous studies or we can construct priors by 
considering no previous knowledge on / so that P(f) = 1 , ex- 
cluding any bias on /. The later will be our assumption for a 
direct comparison between surveys. 

As for any estimation of a random variable, here /, the con- 
fidence interval [/min,/max] in which the true / is can be deter- 
mined with the equation, considering a confidence level CL : 



CL 

If. 



^^"'"^^ P{J\{dj])df (A.3) 



which can be split into implicit equations on /min and /^a, 

1 - CL r^" '" r' 
P(f\{dj})df^ 



^ P(f\{dj})df ^ j P(f\[dj})df (A.4) 



In case of a null detection, Poisson statistic dictates the prob- 
ability of detecting a giant planet around a given star such that 
the likelihood equation lA. 1 [ becomes : 

N 

Null detection sets /min = and the equation I A. 31 becomes 
an explicit equation for given CL : 

-ln{\ - CL) 
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